Hamiltonian structure of thermodynamics with gauge 



Roger Balian 
CEA/Saclay 
Service de Physique Theorique, 
Orme des Merisiers, Bat. 774 
F-91191 Gif sur Yvette Cedex 

and 

Patrick Valentin 
Centre de recherches ELF 
Chemin du Canal, BP 22, F-69360 Solaize 



Abstract 

Denoting by q 1 (i — 1, n) the set of extensive variables which characterize 
the state of a thermodynamic system, we write the associated intensive vari- 
ables 7j, the partial derivatives of the entropy S = S (q 1 , ...,q n ) = qa, in the 
form r y i = —pi/po where po behaves as a gauge factor. When regarded as in- 
dependent, the variables q l ,Pi {i — 0, ...,n) define a space T having a canonical 
symplectic structure where they appear as conjugate. A thermodynamic system 
is represented by a n + 1-dimensional gauge-invariant Lagrangean submanifold 
M of T. Any thermodynamic process, even dissipative, taking place on M is 
represented by a Hamiltonian trajectory in T, governed by a Hamiltonian func- 
tion which is zero on M. A mapping between the equations of state of different 
systems is likewise represented by a canonical transformation in T. Moreover 
a Riemannian metric arises naturally from statistical mechanics for any ther- 
modynamic system, with the differentials dq l as contravariant components of 
an infinitesimal shift and the dpiS as covariant ones. Illustrative examples are 
given. 
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1 Introduction and outline 



In Callen's formulation of thermodynamics jlj, the state of a system is charac- 
terized at each time by the values of n independent quantities, usually extensive 
and conservative, that we shall denote as q l (i — 1, n) . For a single fluid, these 
variables are the internal energy q 1 = U, the volume q 2 = V and the numbers 
q k = N k (k = 3, ...,n) of molecules of each type k; for a pure fluid, we have 
n = 3. For a composite system in local equilibrium, the index i is a compound 
index which refers both to each such physical quantity and to each subsystem; 
for continuous media, these subsystems are volume elements with sizes larger 
than the mean free path. All the equilibrium properties of a given system are 
embedded in the expression of its entropy S = S (a 1 , q n ) , an additive, exten- 
sive and concave function of the state variables q 1 . In particular, the n intensive 
variables 7, are defined as the partial derivatives 

7i=| ? S(/,...,g«). (1.1) 

For a single fluid, they are expressed in terms of the temperature T, the pressure 
P and the chemical potentials per particle fi k as 

7i = ^, 7 2 = 1u = ~y (fc = 3 ' -' n) - (L2) 

The relations (1.1) between the conjugate, extensive and intensive, variables 
constitute the full set of equations of state. 

In non-equilibrium situations, the variables r y i are the local intensive vari- 
ables. The fluxes which govern the dynamics in sufficiently slow regimes are 
expressed in terms of differences of variables 7, for neighbouring subsystems. 
The simplest example is heat transfer across a barrier separating two uniform 
thermal baths. In this case it is sufficient to introduce n — 2 variables, the 
energies q 1 = U 1 and q 2 = U 2 of the two baths, together with the associated 
temperatures r ) 1 = 1/T\ and 7 2 = I/T2. The dynamics is governed by the 
conservation laws 

and by the expression of the heat flux $ across the barrier in terms of the 
response coefficient £(71,72) which characterizes the heat transport through 
this barrier, 

* = -&(7i,7 2 )(7a-7i)- (1-4) 

This formulation of thermodynamics is the one which arises naturally from 
statistical mechanics |l], 0, where the entropy S is identified with the missing 
information, or the disorder at microscopic scale. Equilibrium is then character- 
ized by looking for the largest disorder, subject to constraints on the quantities 
q % . The variables r y i are the Lagrangean multipliers involved in this search. For 
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non-equilibrium thermodynamics microscopic foundations are also provided by 
the maximum entropy method which yields the projection method in the limit 
of negligible retardation effects (see, for instance, the review articles || and 0)). 

Although each state is characterized by the n variables q 1 only, it is con- 
venient for practical purposes and for applications to make use of additional 
variables. One is thus lead to introduce the 2n-dimensional space T with co- 
ordinates q , q n , r y 1 , j n . A thermodynamic system is represented in this 
space by a n-dimensional manifold M characterized by the set of n equations 
of state (1.1); a state of this system is represented by a point in M. For a ther- 
modynamic process, the dynamical equations such as (1.1), (1.3), (1.4) couple 
the 2n time-dependent variables q l and It is here again natural to formulate 

non-equilibrium thermodynamics in the space T, where a dynamical process is 
represented by a trajectory which is constrained by the equations of state (1.1) 
to lie in the manifold M. 

Moreover, these equations of state are often written in an alternative form 
where the entropy and the energy are interchanged. For a single fluid, the func- 
tion S = S(U, V, N k ) is thus inverted into U = U{S, V, N k ). The quantity S can 
therefore be regarded as an additional extensive variable, that we shall denote 
as q° so as to set all the variables q° , , ■■■,q n on the same footing. For a given 
physical system these n + 1 variables are related to one another by 

q Q = S(q\...,q n ). (1.5) 

The dissipation rate is the time-derivative dq°/dt. Altogether, the thermody- 
namic variables q°, q , q n , 7 l7 7„, now appear as the coordinates of a point 

in a 2n + l-dimcnsional space T. In this space, which has a natural contact struc- 
ture [5-8], the dynamical equations involve the n + 1 constraints (1.1), (1-5). 
We have singled out the entropy rather than the energy among the extensive 
variables because our results have a simpler theoretical interpretation if we work 
in the entropy representation (1.5). For practical purposes and for applications, 
it would be easy to transpose the formalism into the energy representation, 
where q° and q 1 are interchanged and where the associated intensive variables 
are T, — P, [i k instead of (1.2). 

The consideration of the space T instead of the initial n-dimensional configu- 
ration space q , ■ ■■,q n has the following interest in equilibrium thermodynamics. 

A physical system is represented in this space T by a rt-dimensional manifold M, 
characterized by the n + 1 relations (1.1), (1.5) or by any other equivalent set of 

n + 1 equations of state. (Strictly speaking M should be called a submanifold 

since it is not defined intrinsically, but as a subset of T.) However experiments 
usually give only indirect indications on the fundamental function S (a 1 , q n ) , 

and are not sufficient to fully determine the thermodynamic manifold M. Since 
less than the n+1 required equations of state are known, the remaining ones 
should be determined by some phenomenological procedure. To this aim, it has 
been advocated |^| to proceed by comparing the physical system under study 
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with some known system having the same degrees of freedom but characterized 
by a different entropy function, or equivalently, by a different thermodynamic 

potential. In the 2n + f-dimensional thermodynamic space T this system of ref- 
erence is described by a known manifold Mo. The incompletely known manifold 
M to be determined is then deduced from M by setting up a correspondence 
between the two systems, represented by some mapping ip in T. The missing 

equations of state for M follow from the corresponding ones for Mo. This map- 
ping can be constructed by means of a continuous deformation scheme. It is 

then generated through infinitesimal transformations in T. We shall regard the 
resulting motion as mock dynamics, and shall interpret the deformation param- 
eter as a fictitious time r. The method can be used to extrapolate, for instance, 
thermal properties from the knowledge of the P, V, T equation of state. In such 
problems we are interested in the structure of a an imaginary dynamical flow 

which would lead from Mo to M, while the physical flow of non-equilibrium 

thermodynamics considered above maps M onto itself. 

The theoretical study of these mappings has already been achieved [8-10]. It 

is based on the remark that the physical manifolds M are not arbitrary. Although 
the relation (1.5) changes from one manifold to another, there exists for each 
one some generating function S which relates the conjugate variables q % and 
7, (i = 1, n) according to (1.1). This is expressed in the 2n + 1-dimensional 

space T by introducing the 1-form 

n 

Co = dq° (1.6) 

i=l 

which is non-degenerate at any point of T. Thus T is endowed with a contact 

structure (T, £j) [5-8]. The existence of some function S (q 1 , q n ) such that 
all the equations (1.1), (1.5) are satisfied is then equivalent to the vanishing of 

(1.6). Among all the n-dimensional manifolds embedded in T, the physical ones 

M must fulfil the condition 



u = (1.7) 

for any infinitesimal variation of a state over M, an exterior differential equa- 
tion which defines Legendre submanifolds. We shall call them thermodynamic 
manifolds to recall that the vanishing of u ensures the fulfilment of the ther- 
modynamic identities (1.1), (1-5). Further conditions on M are imposed by the 
extensivity and the concavity of entropy. 

The transformations M'= ip ^M^j of interest, which map in T a physical 
manifold either onto itself or onto another, should preserve the condition (1.7) on 
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the considered manifolds. More generally, the transformations ip are required to 
map any thermodynamic manifold onto another thermodynamic manifold, that 
is, to preserve the condition (1.7) everywhere. This implies that they are contact 

transformations in T, which multiply the 1-form (1.6) by some non- vanishing 
function A of the coordinates q°,q z ,'~f i . More precisely the push- forward mapping 
ip^ = p^ 1 * which acts on the forms should satisfy 

<p t (u>) = Xu>. (1.8) 

Applications of such contact transformations to thermodynamics have been wor- 
ked out [8-11]. In particular, by means of a continuous set of contact trans- 
formations depending on a single deformation parameter r interpreted as a 
fictitious time, one can explore the equilibrium thermodynamic properties of a 
set of physically different systems labelled by r. 

The structure of the resulting equations of motion is somewhat similar with 
that of Hamiltonian dynamics, though more complicated. However, it differs 
from a symplectic Poisson structure, so that the extension of the standard 
techniques of canonical Hamiltonian dynamics to the present problem is not 
straightforward. The power of such techniques, which for instance readily pro- 
vide variational approaches, makes it desirable to modify the formulation of 
thermodynamics so as to cast the equations of motion into a usual Hamiltonian 
form. We would then benefit both from the flexibility brought in by the idea 
of regarding the 2n + 1 variables q°, q 1 , 7j as independent, and from the whole 
machinery of analytical mechanics. 

We show below that this is feasible at little cost. We first introduce (Section 
2) a non- vanishing gauge variable po, without physical relevance, which mul- 
tiplies all the intensive variables. A new set of variables pi is thus defined as 

Pi = — po7i f° r i = 1) •••,n. The 2n + 1-dimcnsional space T is thereby extended 
into a 2n+2-dimcnsional thermodynamic space T spanned by the variables q t 1 Pi 
with i = 0, 1, n. We associate with a physical system a n+l-dimensional man- 
ifold M in T, parametrized for instance by the coordinates q 1 , q n and p n . A 
gauge transformation which changes the extra variable po while keeping the ra- 
tios Pi/po = — Ji invariant is not observable, so that a state of the system is 
represented by any point of a one-dimensional ray lying in M, along which the 
physical variables q°, q 1 , q n , 7 l5 7„ are fixed. 

We wish to study the transformations in the extended thermodynamic space 
T which map the thermodynamic manifolds M cither onto themselves (in non- 
equilibrium thermodynamics) or onto one another (in the problem of generating 
equations of state). We show (Section 3) that, within a suitable but irrele- 
vant choice of gauge, these transformations are nothing but canonical trans- 
formations of mechanics having some specific features. While the variables 
q l ,"f i (i = 1, ...,n) are conjugate with respect to the entropy in the Legendre 
sense (1.1) of thermodynamics, the variables q l ,Pi (i = 0, 1, ...,n) moreover ap- 
pear as canonically conjugate in the sense of Hamiltonian dynamics. A symplec- 
tic structure is thus induced in the space T. An infinitesimal mapping among the 
manifolds M that preserves the thermodynamic identities is represented by a 
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(possibility time-dependent) Hamiltonian, which generates a flow in T in terms 
of either the real or the fictitious time, depending on the problem. The contact 

transformations in T which relate the physical thermodynamic manifolds M to 
one another are recovered by elimination of the gauge variable pq. 

This relation between contact and canonical transformations is found here as 
a direct outcome of the gauge invariance that we introduced in thermodynamics. 
Actually mathematicians have recognized long ago |l2j that a contact structure 
of the type (1.6) in 2n + 1 dimensions can be embedded into a symplectic 
structure in 2n + 2 dimensions by means of the adjunction of an extra variable. 
This procedure, clearly presented by Caratheodory ||, and later on termed 
as symplectization ]l3fl , has a geometric interpretation in the theory of fibre 
bundles. 

On the other hand, exterior calculus, a standard technique in thermodynam- 
ics Q , has suggested to introduce in the 2n-dimensional space of extensive and 
intensive variables a symplectic structure which establishes a duality between 
these variables. For instance, in the energy representation |hJ, one currently 
considers for a fluid the symplectic 2-form dT A dS — dP A dV in the corre- 
sponding 4-dimensional space. Likewise, in the entropy representation pH] , the 
fundamental symplectic 2-form is d (1/T) A dU + d (P/T) A dV. More generally, 
the 2n-dimensional thermodynamics space T can be endowed with a symplectic 

n 

structure generated by the 2-form A dq 1 . The vanishing of this form is 

i=l 

known to characterize the surface M of T which describes the equations of state 
(1.1) of any given thermodynamic system, so that M is a Lagrangean submani- 
fold of T. Deformations which map thermodynamic systems onto one another in 



the 2n-dimensional space T have also be considered 16 . They will be recovered 



below (see Section 4) from canonical transformations in our 2n + 2-dimensional 

n 

space T endowed with the symplectic form ^2 dpi A dq 1 , as special cases when 

i=0 

the entropy variable q° may be left aside and when the gauge may be fixed as 
Pq = — 1. Our Hamiltonian dynamics in T thus generalizes and encompasses 

both existing formulations in the 2n- and 2n + 1-dimensional spaces T and T. 

The present formalism applies in particular to the dynamical equations in 
non-equilibrium thermodynamics (Section 4), which can thus be cast into a 
Hamiltonian form. Given a set of equations of motion, it is usually not obvious 
to recognize whether they have a Hamiltonian nature. This was done for the 
hydrodynamics of ideal fluids Jl7t . Here we find another type of Hamiltonian 
structure, for any dissipative system in the thermodynamic regime. The exis- 
tence of this structure is based on the idea that any dynamical system can be 
embedded into a Hamiltonian system with a double dimension (see ref. Jl5| , 
Chap. 10). The remarkable fact in non-equilibrium thermodynamics is the ex- 
istence of a direct physical interpretation for the additional variables. In case 
the entropy is left aside, the evolution of the n independent variables q , ■■■,q n 
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which characterize the state of the system at each time has the same struc- 
ture as in analytical mechanics in spite of the existence of dissipation. If wc 
regard the set q 1 , q n as position variables, their associated intensive variables 
7 1; ...,7„ appear as their conjugate momenta, and the dynamics is governed by 

an effective Hamiltonian / in the 2n-dimensional space T. In case the entropy in 
included among the dynamical variables, a similar canonical structure is found 
provided we pay the price of introducing the gauge variable po upon which the 
effective Hamiltonian h should depend. Whereas the conjugate momentum of 
each q % is now pi for i = 1, ...,n, the conjugate momentum of the entropy q° 
is po, and the dissipation dS/dt is expressed as dh/dpo. The constraints (1.1), 
(1.5) satisfied at the initial time are conserved by the Hamiltonian flow, which 
moreover implies that / — or h — along the physical trajectories. The fact 
that only the part of the Hamiltonian flow located in the manifold (1.1), (1.5) 
is relevant for physics solves the paradox of the existence of a Hamiltonian for 
dissipative dynamics. As an illustration we work out in Section 4 Hamiltonian 
formulations for the equations of motion (1.3), (1.4) describing heat transport. 
We show moreover that these results are readily extended to any dissipative 
process in the thermodynamic regime, whether mechanical, thermal, electrical 
or chemical. 

The equations of state for different physical systems are also generated from 
one another by means of Hamiltonian transformations in the space T. We illus- 
trate this by writing in Section 5 Hamiltonian mappings between different van 
der Waals fluids. 

Finally we recall in Section 6 how quantum statistical mechanics generates 
[|| a natural metric structure in thermodynamics |l^, [l^] . The thermodynamic 

Lagrangean manifolds M, M or Legendrean manifolds M thus also acquire a 
structure of Riemannian manifolds. The metric is obtained as 



and the thermodynamically conjugate variables now appear as contravariant and 
covariant coordinates on the manifold. The concavity of the entropy function 
associated with M is equivalent to the positivity of ds 2 over M. This property 
may be used to select the physically admissible manifolds, which should not 
only satisfy the algebraic condition (1.7) but should also have a positive ds 2 so 
as to describe thermodynamically stable systems. 

2 Thermodynamics as a gauge theory 

We start with a remark, drawn from statistical mechanics, about the definition 
of the intensive variables. Consider, for instance, a fluid in grand canonical 
equilibrium. In terms of the intensive variables r y i defined by (1.2), its density 




(1.9) 
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operator D (in quantum statistical mechanics) is expressed by 



D oc exp 



1 / 



(2.1) 



where H and N k are the Hamiltonian and particle number operators. The oc- 
currence of Boltzmann's constant fce arises from the choice of the kelvin as the 
unit of temperature and the joule per kelvin as the unit of entropy. We note 
from (1.2) that all the 7,'s are inversely proportional to the temperature. Thus, 
multiplying both fee and the 7^ ' s by a constant does not affect the physics, since 
it amounts to changing the unit of temperature, or equivalently the irrelevant 
coefficient in the definition of the von Neumann entropy in quantum statisti- 
cal mechanics, which is identified at equilibrium with the absolute entropy of 
thermodynamics within a multiplicative constant. 

This suggests us, for an arbitrary thermodynamic system, to introduce an 
additional variable po and to replace the set j t by the new, intensive, scaled 
variables 



Pi = -poli (i = l,...,n). (2.2) 

In the resulting 2n+2-dimensional space T, a dilation of the variables q l , keeping 
the variables pi unchanged, is a physical operation which leaves the equations 
of state (1.1), (1-5) unchanged in case the system is extensive. Symmetrically 
we introduce a mathematical operation, the dilation of the intensive variables, 
keeping the extensive ones unchanged, 

Pi^Xpi, q i ^q i (i = 0,l,...,n), (2.3) 

where A is a non-zero constant. This operation, which does not affect the phys- 
ical variables q , ...,q n ,j i , ...,7„, can be regarded as a gauge transformation of 
the first kind. In the trivial gauge po = —1, we have pt = 7$. A change of gauge 
which leads to po = —1/ks suppresses Boltzmann's constant from the density 
operator (2.1) while changing the 7^'s into PiS. (The macroscopic entropy re- 
mains, however, unchanged as its expression in terms of D still contains a factor 

We regard below, more generally, the dummy factor po in (2.2) not as a 
constant but as an independent variable. We thus allow gauge transformations 
(2.3) of the second kind, where A is some function of the 2n + 2 coordinates 
q l ,Pi- An illustration will be given by Eq. (2.10), which is obtained from the 
trivial gauge po = — 1 by taking A = po/pi- 

In this formulation of thermodynamics as a gauge theory, the variables po,Pi, 
...,Pn arc not completely meaningful physically, whereas their ratios r y i given by 
(2.2) remain invariant in a gauge transformation. The situation looks like clas- 
sical electromagnetism, where the e.m. potential is modified in a gauge trans- 
formation and is thus physically unobservable (as the variables Pi), whereas the 
e.m. field is well-defined (as the physical intensive variables 7J. We shall work 
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in the 2n + 2-dimcnsional extended thermodynamic space T with coordinates 
q\Pi, including the gauge factor pq. This is mathematically convenient, and 
physical interpretation will require to return to the 2n + 1 physical coordinates, 
namely the entropy q°, the n extensive variables q % and the n intensive variables 
7 i of Callcn. 

In the extended thermodynamic space T, a state for a system is not repre- 
sented by a single point, but by any point of a one-dimensional ray characterized 
by fixed values of the coordinates q°,q 1 ...,q n and of the ratios r y i = —pi/po- The 
points along this line, equivalent for physics, result from one another by a gauge 
transformation (2.3). A given system is described by a n + 1-dimcnsional sub- 
manifold M, which may be parametrized by the n extensive variables q , q n 
and by the gauge fixing variable po- The entropy function S(q 1 , ...,q n ) of the 
system generates the remaining n + 1 variables as 



q = S(q\...,q"), P, = -P0 dS{ql dq r qn) - (2.4) 

We recalled in the introduction that the existence of such a generating function 

S characterizes the thermodynamic submanifolds M in the 2n + 1-dimensional 

space T. It is then expressed by the vanishing (1.7) of the 1-form (1.6). Accord- 
ingly, in the extended 2n + 2-dimensional space T, we introduce the 1-form 



n 

uj = ^2p l dq\ (2.5) 

i=0 

and the n + 1-dimensional thermodynamic manifolds M are characterized by the 
vanishing of this form: 



u = over M. (2.6) 
The 1-form u> induces a symplectic structure 

n 

duj = ^2 d Pi A d( i ( 2 - 7 ) 

i=0 

on T that we denote as (T, dui) . Any thermodynamic manifold M belongs to the 
set of the so-called Lagrangean manifolds in T, which are the integral subman- 
ifolds of du> with maximum dimension (n + 1) . Moreover M is gauge invariant, 
which is implied by (2.6) but not by dui = 0. 

When the set q % contains all possible extensive variables, the extensivity of 
the entropy function S(q 1 , q n ) is expressed by the Gibbs-Duhem relation 

i—1 
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Translating it into the space T by means of Eq.(2.4), we find 

n 

X> = °- ( 2 - 9 ) 

i=0 

This condition defines the 2n + 1-dimcnsional extensivity sheet E in the space T, 
and the thermodynamic manifolds M should lie in this case within this surface 
E. However, if the systems are restricted to contain a fixed quantity of matter 
or to have a fixed volume, the corresponding manifolds M are free from the 
above constraint. 

Statistical mechanics exhibits the different nature of the entropy q° and of the 
other extensive variables q , ■■■,q n . However it is traditional in thermodynamics 
to introduce a symmetry between all q l, s for i = 0, 1, ...,n. In particular the 
energy as function of the entropy, rather than the converse, is often used as 
a thermodynamic potential. This symmetry is conveniently reflected in the 
present formalism by the use of the gauged intensive variables po,pi, ...,p n , as 
shown in (2.5), (2.6) and (2.9). For instance, for a fluid, with the choice of gauge 
Po = — 1, the intensive variables pi, ...,p n reduce to the intensive variables j i 
arising from the entropy representation, given by Eq. (1.2). However, with the 
choice of gauge p\ = —1, the other intensive variables pi, namely 

Po = — = T, p 2 = -pol 2 = - p , Pk = -Polk = ( k = 3, -,n) , (2.10) 

are those which arise from the energy representation of thermodynamics, where 
the relation q 1 = U {q° , q 2 , ...,q n ) replaces (1.5). Thus, the gauge transformation 
(2.3) with A = po/pi = — 1/7 X amounts to a switch from the entropy to the 
energy representation. 

3 Canonical transformations in the extended 
configuration space 

We proceed to study the transformations which, in the extended thermody- 
namic space T, map the thermodynamic manifolds M onto one another. Such a 
mapping may describe a non-equilibrium thermodynamic process, in which case 
the initial and the final manifolds are the same, or it may describe a mapping 
between the equations of state of different physical systems. 

The thermodynamic submanifolds M satisfy the condition (2.6). We impose 
that this condition should be preserved by the considered mappings M' = ip(M) 

for any M. We recalled in Eq. (1.8) that, in the 2n + 1-dimensional space T 

the corresponding mappings M' = (p(M) should satisfy <I>' = <^„(w) = Xui where 
A is some non- vanishing function of the coordinates. Likewise the mappings 
M' = ip(M) in the 2n + 2-dimensional space T that preserve the thermodynamic 
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identities should satisfy 



ui' = <p*(u>) = \ui, (3.1) 

a property which ensures that us = implies us' = 0. However, the factor A in 
(3.1) can be absorbed by a gauge transformation, since the multiplication of us 
by a function A of the coordinates q l ,Pi can be achieved by the gauge transfor- 
mation (2.3) involving the same factor A. Hence, within an irrelevant change of 
gauge in M', any mapping M' = ip(M.) which satisfies (3.1) is equivalent to a 
transformation in T which conserves the 1-form us. We can therefore take advan- 
tage of the gauge invariance and restrict ourselves to the mappings M' = ip(M) 
such that 

us' = tp*(u}) = us. (3.2) 

The condition (3.2) implies that the 2-form (2.7) obtained by taking the 
exterior derivative of (2.5) is also conserved by the mapping tp, that is, 

ip if (duj) — dus. (3-3) 

This means that, if we endow the space T with a symplectic structure with 
Poisson brackets 

{'/'•/',,} S), (3-4) 

the considered mappings are canonical transformations of the same type as in 
analytical mechanics ^ [l3| p0[ . 

Not every canonical transformation is admissible, however, since the conser- 
vation (3.2) of the 1-form us is a stronger condition that the conservation (3.3) 
of the 2-form dus. Consider an infinitesimal canonical transformation, generated 
by the (possibly time-dependent) Hamiltonian h(q°,q l , q n ,po,pi, ■ ■■,p n ) '■ 

dh . dh 

The evolution takes place either as function of the physical time t in non-equili- 
brium thermodynamics, or as function of a fictitious time r for continuous defor- 
mations relating different systems to one another. The conservation of us reads 



. , . , dh dh , 

p t dq +pidq I = ) p, t d- ^-rdq 



i=0 i=0 



dp t 




(3.6) 



for any variation of the coordinates q' ,pi (i = 0, 1, n). As h is defined within 
an additive constant, we can impose that it satisfies the identity 
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or equivalently 

h (q ^ 1 , ...,q n , \po, Xpi, -,Xp n ) = A/i (g°,<7\ ...,q n ,p ,pi, -,p n ) , (3-8) 

where A is a constant or more generally a non- vanishing function. The condition 
(3.8) also ensures that, if two points are deduced from each other by a gauge 
transformation (2.3), their images in the evolution remain related by the same 
transformation. The rays which represent the states in the space T are thus 
deduced from one another in the Hamiltonian motion, which therefore preserves 
the gauge invariance of the manifolds M in the dilation (2.3). Actually we 
have noted that the thermodynamic manifolds M (for which u> = 0) are n + 
1-dimensional Lagrangean submanifolds of T (for which dtu = 0) which are 
moreover gauge invariant. 

Altogether the mappings in the extended thermodynamic space T which 
preserve these two properties that characterize the thermodynamic structure of 
the manifolds M are canonical transformations of analytical mechanics, gener- 
ated by a Hamiltonian h which is a homogeneous function with degree 1 in the 
variables pi . 

If the set q°, q , q n includes all the extensive variables, the mappings 
should leave the extensivity sheet E defined by (2.9) invariant. Expressing the 
time-derivative of (2.9) by means of (3.5) and using (3.7), we obtain in this case 

i=0 H 

Hence, h should also be a homogeneous function with degree 1 of the variables 
q % . Note that, in spite of their formal analogy, the identities (3.7) and (3.9) have 
different status, since dilation of the p^s is unphysical while dilation of the q l, s 
is associated with extensivity, and since (3.7) should always be satisified. 

A finite canonical transformation may be characterized by its generating 
function H(q , q n ,P(,, ...,p' n ), which produces the mapping q\pi i — ► q h ' ,p\ as 

Pi = ^, <f=™ (i = 0,1, ...,»). (3.10) 

(If the Jacobian of this mapping vanishes, one should use other generating func- 
tions, where some initial variables q % are replaced by PiS, and conversely for the 
final variables; see ref. (l3); an example will be given by Eq. (5.8) below.) The 
conservation of the 1-form u> is expressed, through the same calculation as in 
(3.6), by 

so that H is a homogeneous function of the variables p\ with degree 1 . Infinites- 
imal transformations (3.5) are recovered from 

n 

H*J2Pi <i + £h{q ,...,q n ,p' Q ,...J n ) (3.12) 
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for e small. Note that the condition (3.11) results in the vanishing of a Legendre 
transform of TC with respect to the whole set of variables p\. Hence no generating 
function depending on the variables q l ,q n can produce here the Hamiltonian 
mapping. Extensivity implies, as in (3.9), that TL is a homogenous function of 
the variables q l with degree 1. In this case, the identity (3.11) implies that 
H vanishes when the points {q l ,Pi) and (q n ,p' i ) which are in correspondence 
through (3.10) lie on the extensivity sheet. 

We recover contact transformations in the space T by the projection n which 
eliminates the gauge factor pp from Hamilton's equations, a procedure inverse 
from symplectization || Indeed, the time-dependence of any function 

g [q°, q n , r y 1 , 7 n ) , where the 7j's are expressed in terms of the canoni- 
cal variables by (2.2), is given, according to Hamilton's equations (3.5), by the 
usual Poisson bracket 

9 = i9M = ^ WWi - WWi , (3-13) 
generated by (3.4) in the space T. We can calculate it by using the fact that 

f(q°, <f , 7l , ..., 7 „) = --h (q°, ...,q n , Po , ..., Pn ) (3.14) 

is a function of the variables 7, ; = —pi/po only, as a consequence of the homo- 
geneity condition (3.8). The result, 

• d 9 ( , d f\ sr^ dg df ^dgfdf df\ , oir . 

9 = 2& {-f + E^J + E - E Wi [ Wi + 7^ J - (3-15) 

generates the known contact flow || in the thermodynamic space T. The greater 
simplicity of the canonical flow in the extended space T arises from the mock 
dynamics 

Po = (3.16) 

that we have introduced for the gauge variable, which allowed us to get rid 
of the factor A in (3.1), whereas contact transformations should manage with 

this factor in (1.8). We can thus regard a contact transformation in T as a 
Hamiltonian transformation in T followed by the gauge transformation which 
brings po to its initial value —1. However, the latter gauge transformation (2.3) 
cannot in general be generated by a Poisson structure, so that this structure is 
lost in the elimination of pg. As it often occurs in gauge theories, the equations 
become more complicated when the gauge is fixed without care. 

Note that the antisymmetry of the Poisson bracket (3.13), which gener- 
ates the motion in the 2n + 2-dimensional space T, is lost when this motion 

is projected as (3.15) on the 2n + 1-dimensional space T. Indeed the first term 
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— (dg/dqo) / of (3.15) has no counterpart. We can trace back this lack of symme- 
try to a hidden difference which already existed between g and h in the Poisson 
bracket (3.13). Being a physical quantity, g is a homogeneous function of the 
variables pi with degree 0. However, h, the generator of the motion in the space 
T, is not a function of the physical variables q l , j i only, but it includes po and is 
homogeneous in the p^s with degree 1. The first term of (3.15) arises from this 
difference of behaviour of g and h in a gauge transformation. 



We can represent the relations between the contact space (T, Q) and the 
symplectic space (T, dw) by the following diagram: 



The symplectization, usually performed fliaj by constructing T as a fibre bundle 
associated with the configuration space cp, q n , amounts here to the intro- 
duction of a gauge structure in the space of intensive variables j i = —pi/po, 
whereas the projection tt amounts to the fixation of the gauge through po = — 1. 
The Hamiltonians h which are equivalent to contact transformations / are con- 
strained by the homogeneity condition (3.8), which ensures that the gauge in- 
variance of the Lagrangean manifolds M is preserved. 

4 Hamiltonian equations for dissipative dynam- 
ics and quasi-static processes 

The above formalism holds in particular for a continuous mapping which leaves 
some physical manifold Mo invariant. The equations of motion (3.5) can thus 
describe non-equilibrium thermodynamic processes which take place in a given 
system characterized by Mo, and which are represented by the motion of a 
point in Mo. In these dynamics, the other thermodynamic manifolds M are 
not kept invariant but are transformed into one another, contrary to Mo- The 
effective Hamiltonian h will therefore depend on the equations of state (2.4) 
which parametrize Mo (or equivalcntly on the entropy function S of the system), 
as well as on the transport coefficients. This effective Hamiltonian h should not 
be confused with the microscopic Hamiltonian H which governs the dynamics 
in statistical mechanics. 

Apart from its general property (3.8), h should here be such that the Hamil- 
tonian flow that it generates through (3.5) lies, for any point of Mo, in Mo 
itself. In particular, the relation q° = Stq 1 , ■■■,q n ) should be conserved, which 
is expressed by writing its time-derivative as 



(T, duj) 
Symplectization j 



symplectic transformation, h 



(T, dtu) 

I Projection, tt 



contact transformation, / 




i=l 



i=i 



(4.1) 
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Together with (3.7), this implies that 



h = overM . (4.2) 

Conversely, the general homogeneity property (3.8) of h and the vanishing of 
h on the manifold Mo which represents the given system in the extended ther- 
modynamic space T are sufficient to ensure that the equations of state (2.4) 
remain satisfied at all times provided they are satisfied at the initial time. This 
is readily checked by regarding the set (2.4) as a parametrization of M in terms 
of the independent variables q , ...,q n ,po- This can also be seen as a consequence 
of (4.2) and of the conservation, along the motion, of the 1-form u> which thus 
remains zero. The vanishing of h over Mo is the counterpart of the vanishing of 

/ in the corresponding contact transformation in T that keeps Mo invariant, a 
property rigorously proven in ref . ^] . 
One among the equations of motion, 

i° = JT> (4-3) 
dpo 

describes the rate of change of the entropy, which for an isolated system is the 
dissipation. Thus, although the variable po has in itself no physical meaning, 
its occurrence in the effective Hamiltonian is essential to deal with dissipative 
processes. It would clearly not have been possible to assign to an irreversible 
process Hamiltonian equations including the time-derivative q of entropy with- 
out introducing such an extra variable. Dynamical features which apparently 
contradict usual properties of Hamiltonian motions, such as the fact that the 
dissipation (4.3) cannot be negative, or the convergence at large times of the 
trajectories towards a fixed point which describes global equilibrium of the sys- 
tem, arise from a special choice of the initial point, which should always lie 
on Mo. In particular, the Liouville theorem is compatible with the convergence 
of physical trajectories towards a fixed point, because there exist neighbouring 
unphysical trajectories, close to Mo but outside it, which diverge away from this 
fixed point. We shall illustrate this fact in Eq. (4.13) below. 

Since the Hamiltonian equations (3.5) are here physically meaningful only 
for the flow in M , the Hamiltonian h is not defined in a unique fashion for a 
given physical process. Changes of h which do not modify its value (h = 0) on 
M and its first-order derivatives on M are irrelevant for physics. 

As an illustrative example, let us consider heat transfer between two thermal 
baths. We assume the thermal conductivity of each bath to be much larger 
than that of the barrier, so that their two energies U 1 and U 2 are sufficient 
to characterize at each time the state of the system. The entropy S is the 
sum of the entropies Si (U 1 ) and 5*2 (U 2 ) of the two baths. We recalled in 
the introduction the equations of motion (1.3), (1.4) for this model. In the 6- 
dimensional extended thermodynamic space T, the manifold Mo, parametrized 
by q 1 — U 1 , q 2 — U 2 and the gauge variable po, is characterized by the equations 



'1 



= S 1 (q 1 ) + S 2 (q 2 ), p^-poP^q 1 ), p 2 = - P oP 2 (q 2 ), (4.4) 
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where we denote as 



M u * = ^mr> fi *™ = ^mr (4 - 5) 

the inverse temperatures of the baths regarded as functions of their energies. 
We readily check that on the manifold (4.4) the equations of motion (1.3), (1.4), 
rewritten in terms of q l ,Pi (i = 0, 1, ...,n), can be generated through (3.5) from 
the Hamiltonian 



(4.6) 



which depends on q 1 and q 2 through (3 1 (q 1 ) and /3 2 (q 2 ). As it should, h vanishes 
on the physical manifold (4.4), and is of degree 1 in the variables pi. The latter 

n 

property ensures that the dissipation (4.3) is the one, 7i9*i associated with 

i=l 

the flux (1.4). The conservation of energy is reflected by the fact that h depends 
on pi and pi only through their difference. 

As noted above, we can alternatively use for h any Hamiltonian equivalent to 
(4.6) to first-order in the variables p\ +po f3 1 (q 1 ),p2+Pofl2(q 2 ) aim Q° ~ ^(q 1 7 q 2 ) 
which vanish on the manifold Mo. For instance, we can take, in terms of the 
flux $(7i,7 2 ) = L(j l7 -f 2 ) (7 2 — 7 X ) of energy escaping from the bath 1, which 
is a function of the intensive variables, the Hamiltonian 

h = *(--, --) [(pi + P0 (3 1 ) - (pa + p (3 2 )] ■ (4.7) 
Po Po 

The Hamiltonians (4.6) or (4.7) do not depend on the entropy variable q°, 
so that the time-derivative 

dh , , „. 

P» = - W (4-8) 

of the gauge variable po vanishes. The elimination of the gauge is trivial in 
such a situation. By fixing the gauge so that po = —1 at all times, and by 
keeping aside the entropy variable q°, we can then regard the 2n-dimensional 

space T spanned by q , q n , 7 l5 7„ as a mechanical phase space in which the 
extensive and intensive variables q % and ■y i (i = 1, ...,n) appear as canonically 
conjugate. Their dynamics are generated by the Hamiltonian p f = —f defined 
by (3.14), as 

q l = ^- <y. = -^ (4 9) 

Actually Eq. (3.14) defines in general a contact Hamiltonian / in the 2n + 1- 
dimensional space T, but the absence of q° allows us to regard / as an ordinary 
Hamiltonian in the 2n-dimensional symplectic space T, in the special case when 
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the effective Hamiltonian in the space T does not depend on the entropy variable 
q°. For the heat transfer problem, the Hamiltonian / which arises from (3.14) 
and (4.6) is, in terms of the "position" variables q x ,q 2 and their conjugate 
"momenta" 7i,7 2 , 

/ = \L (fi lt 2 ) [( 7l - 7 2 ) 2 - (A ~ P2f] • (4.10) 
Equivalently as regards the physical dynamics, we get from (4.7) 

/ = * (7i,7 2 ) [(7i - Pi) - (7 2 - &)] • (4.11) 



The Hamiltonians (4.6) in T or (4.10) in T involve two terms, the first one 
of "kinetic" type, the second one of "potential" type. These two terms occur 
with opposite signs, in contrast to what usually happens in mechanics. For 
the motions in which we are interested, they exactly cancel out at the initial 
time, and hence at all times. As indicated above, the special choice of an initial 
point in Mo prevents the divergence of trajectories which is expected from the 
wrong sign of the potential term. We illustrate this point by fully solving the 
equations (4.9) in the whole space T for a model where L is a constant and 
where the entropies are quadratic functions of the energies: 

Si(ir)=a i U i -\b i {U i f , (i = 1,2). (4.12) 

(The Gibbs-Duhem homogeneity condition (2.8) is not satisfied here because 
only one extensive variable, q l = U l , has been introduced for each subsys- 
tem; hence b^ 1 is an extensive quantity.) The Hamiltonian (4.10) is that 
of a two-dimensional harmonic oscillator in which the sign of the potential 

— \L (p 1 — P 2 ) 2 = — \L ( a i — ^l? 1 — a 2 + hq 2 ) has been inverted. Its general 
flow is 

7 1 -7 2 = Ac rt + Bc- rt , 0! -(3 2 = -Ac Tt + Bc- T \ (4.13) 

where F = (b\ + b 2 ) L, with the total energy q 1 +q 2 and the quantity 6 27l + &i 72 
as constants of the motion. The concavity of the entropy functions (4.12) implies 
bi > 0, and the positivity of the response coefficient L then implies T > 0. The 
arbitrary constants A and B are determined by the initial conditions, which 
for a physical process entail 7l — 7 2 = f3 1 — [3 2 and hence A = 0. Thus only 
the decaying terms Be~ rt occur in (4.13) for physical processes, although the 
diverging terms Ae rt are present for non-physical trajectories which lie outside 
the manifold Mo- Likewise the dissipation evaluated from (4.12) and (4.13), 
S = L (i? 2 c~ 2r * — A 2 e 2rt ) , remains positive at all times only on M . It is 
therefore important, in order to use Hamilton's equations (3.5) or (4.9) for a 
numerical solution of dissipative motions, to enforce the constraints (2.4) which 
characterize Mo, although these constraints are conserved along the motion, 
because small errors may produce increasing spurious effects. 

In spite of its simplicity, the heat transfer problem considered above is a pro- 
totype for any process of non-equilibrium thermodynamics, describing transfer 
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of heat, of momentum, of particles, or chemical reactions. In any such case, the 
time-derivative of each extensive variable q % , where the compound index i = k,a 
refers both to the nature k of the quantity transferred and to the considered 
subsystem a, is expressed as a sum of outgoing fluxes (71, •■■,7„) , which 
depend on the intensive variables r y i and which involve the subsystems (3 that 
interact with a. The conservation laws are expressed by = — $^' 3 . If we keep 
aside the entropy variable q° and fix the gauge as po = — 1, we can generate the 

dynamics in the 2n-dimensional space T from the effective Hamiltonian 



f = 52**13 (7l,-,7n) 



dS (g 1 ,-,^) 



7i 



(4.14) 



This expression generalizes (4.11), which involved only two opposite fluxes 
$ x = $ = — $ 2 . It also applies to non-isolated systems, in which case some 
fluxes are imposed by external sources. Since the trajectories of interest satisfy 
the constraints j i = dS/dq l , the Hamiltonian flow in the whole space can be 
stabilized around those constraints by adding to (4.14) terms proportional to 



OS {q\. 



dq % 



(4.15) 



which do not affect the dynamics on the physical manifold M. 

Effective Hamiltonians such as (4.14) are not arbitrary. They should satisfy 
several properties imposed by the theory of non-equilibrium thermodynamics 
[[[} ^|. The entropy function S {q 1 , q n ) should be concave, expressing stability 
of matter. The fluxes should be two by two opposite, so as to ensure the 
conservation laws. The matrix of response coefficients L should be positive so 
as to ensure that the dissipation is not negative. Finally the symmetry and 
invariance laws should be reflected in the form of the functions S and $. 

In the special case of quasi-static processes, there is no dissipation. Such 
processes are usually considered for a non-isolated system which remains nearly 
in equilibrium at each time under the effect of external sources. If it does not 
interact with a heat source, its entropy remains constant, so that the effec- 
tive Hamiltonian h should not depend on pq. Consider, for instance, a fluid 
in adiabatic expansion. Its instantaneous state is characterized, for a fixed 
particle number, by the two variables q 1 = U and q 2 = V, and its equi- 
librium properties by the entropy q° = S (U, V) . The intensive variables are 
7j = —pi/po = 1/T, 7 2 = —P2/P0 = P/T. If the motion of a piston changes the 
volume sufficiently slowly, according to V = the flux of energy is U = — P$ 
and the resulting quasi-static dynamics of the fluid in the 6-dimensional space 
T can be generated by the effective Hamiltonian 



h = $ 



dS IdS 

P1 W mJ- p2 



(4.16) 
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(which possibly depends on time through <&). The constancy of entropy is ob- 
vious from dh/dpo = 0. In the 4-dimensional symplectic space T = U, V, 7 l5 7 2 , 
the resulting effective Hamiltonian is 



a simplified form of (4.14). 

For continuous media, the index i in q % and r y i for i ^ not only refers 
to the nature of the variable but also labels the volume elements. This index 
thus includes the coordinates in ordinary space, so that the variables q % and 
7, constitute fields (such as the energy density or the local velocity). The 
expression (4.14) becomes the effective Hamiltonian for a canonical field theory. 
Moreover, in hydrodynamics, the fluxes need not vanish with the gradients of 
the intensive variables. 

We have noted that, under the condition (4.2), the dynamics generated by 
the Hamiltonian h (or /) preserve the constraints (2.4) which express the equa- 
tions of state, provided they are satisfied at the initial time. One could take 
advantage of these constraints so as to reduce the number of variables, or to 
obtain by standard techniques non-canonical though Hamiltonian dynamics in 
terms of reduced Poisson structures Q. 

5 Generation of equations of state 

We have written the equations of state which determine a thermodynamic man- 
ifold M in the specific form (2.4). In practice a thermodynamical system is 
characterized in the 2n + 2-dimensional space T by generalized equations of 
state, which constitute a set of n + 1 equations equivalent to the set (2.4). As 
a consequence of gauge invariance, each of them is homogeneous in the pi's. 
These equations of state are not always well known, and the missing ones may 
be determined phenomenologically by comparison with known similar systems. 



To this aim, contact transformations in the space T have been used [8-11]. We 
suggest here to rely on canonical mappings or Hamiltonian flows in the space T 
in order to transform the sets of equations of state of different system into one 
another. 

As an illustration, let us consider two systems characterized by their entropy 
functions S (q 1 , q n ) and S' [q 1 , q n ) , respectively. We wish to map their 
associated thermodynamic manifolds M and M' onto each other in T. For 
simplicity, we consider a mapping which does not affect the coordinates q , q n . 
The generating function TL which achieves the canonical mapping through Eqs. 
(3.10) is not defined in a unique fashion, since its action on manifolds other 
than M and M' is not specified. We can readily check that the following choice 




(4.17) 



19 



is suitable: 



H = Y,P^+Po i S ' <?") - S (q\ q n )] ■ (5.1) 



Indeed, Eqs.(3.10) yield 



OS dS' 
Pa = P6, Pt+Po d^ = p 'i +P '°dqT> 



= go -S, < f=q i (i = l,...,n), (5.2) 

which are obviously satisfied for a pair of points located on M and M' and having 
the same coordinates poj q\ ■ q n ■ Eq.(5.1) exhibits the occurrence of a thermo- 
dynamic potential, here the entropy, in the expression of the generating function 

Ti.. Such an occurrence was already recognized for contact transformations in T 
[9-11]. 

For a continuous set of hypothetical systems labelled by a deformation pa- 
rameter r regarded as a mock time, the representative manifolds M(r) in T 
can be parametrized according to (2.4) by means of the family S T (q 1 , ■■■,q n ) 
of entropy functions. The flow which transforms them into one another with- 
out changing the coordinates po, q 1 , q n is generated, according to (3.12) and 
(5.1), by the time-dependent Hamiltonian 

h = P^S T (q\...,q n ). (5.3) 

More general flows transforming continuously the manifolds M (r) into one an- 
other can be generated by Hamiltonians, obtained by adding to (5.3) some 
r-dependent function which vanishes on M (r) and is homogeneous with de- 
gree 1 in the p^s. Indeed, as shown in Section 4, this additional term lets the 
coordinates po^q 1 , change in time without modifying M (t) . 

In the above example, we have assumed the entropy to be given as a ther- 
modynamic potential. For application purposes we can take advantage of the 
geometric nature of the formalism, which allows us to parametrize the thermo- 
dynamic manifolds with variables other than q , q n ,po and accordingly to use 
different thermodynamic potentials. As an example, let us reconsider transfor- 
mations which map a van der Waals fluid onto another one or onto a perfect gas, 
a problem already studied in the framework of contact transformations |ll[ . 
We keep here the particle number fixed: we shall take for it the Avogadro num- 
ber N\. The space T has thus 6 dimensions. The practical variables, namely, 
the molar entropy, energy and volume, the temperature, the pressure and the 
single-particle chemical potential /1 obtained from (2.8) are identified as 

S = q°, U = q\ V = q\ T = — — , P = ^, 

Pi Pi 

flN A = U -TS + PV = — ( Po q Q + P iq 1 + P2 q 2 ) ■ (5.4) 

Pi 
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The occurrence of the gauge factor pi in the denominators of the intensive 
variables instead of po is related to the use of the energy representation, more 
convenient here than the entropy representation. In fact, if we take the free 
energy F (T, V) as a thermodynamic potential instead of S (U, V) , we find the 
generalized equations of state for a manifold M in the form 



q» = S = - 



dF 



q i + Pl q ° = U-TS = F, * 
Pi Pi 



= P = - 



(5.5) 



where the arguments in F and its derivatives are replaced by —po/pi = T and 
q 2 = V. Here, M is thus parametrized in terms of q 2 ,po and p\. 
A van der Waals fluid is characterized by its molar free energy 



F (T, V) = -RT\n (V - 6) - — + * (T) , 



(5.6) 



where R = Aa&b and where the constants a, b and the function 'J (T) depend 
on the fluid. An ideal gas of structureless particles with mass m corresponds to 
the special case 



6 = 0, y(T) = RTln 



Aa ( 2irh 2 \ 3/2 ' 
~~c~ \mk B Tj 



(5.7) 



We wish to map this van der Waals fluid onto another one (or onto an 
ideal gas), characterized by its free energy F' (T, V) of the form (5.6) where 
a, b, (T) are replaced by a', 6', \&' (T) . We consider in T a mapping of M onto 
M' for which the coordinates q 2 ,po and p\, that is, T and V, remain unchanged. 
This is achieved by introducing a generating function K,, the Legendre transform 
of TL with respect to p[ and q 1 , due to the choice of variables T and V in the 
thermodynamic potential. Transposing Eq.(5.1) yields 



ic(q , pi, q 2 , Po- <i n , P' 2 )= Pol -Pil' 1 +P2Q 2 +Pi 



pi J V pi 



(5.8) 



where F and F' have the form (5.6) or (5.7). We readily check that the mapping 
provided by the partial derivatives of JC in the full space T reads: 



q° + 



Po 

dF 
df 



Po- 



Pi =Pl, 

dF' 
~dT 



dF , ,dF' , 2 
P2+Pi W =P 2 +P 1W , q 



q\ 



1 PodF n p> dF> 

- t = q - 



q - 



Pi dT 



P'l dT 



- F'. (5.9) 



Hence the equations of state (5.5) for M imply the corresponding ones for M', 
and /C therefore maps as expected these manifolds onto each other. 

If a, 6, (T) depend continuously on a mock time r, the free energy (5.6) 
depends on r, what we denote by F T . The generating function K, of (5.8) taken 
between the times r and dr yields a time-dependent Hamiltonian, as TL does in 
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(3.12). The family of manifolds M (r) is then obtained from the Hamiltonian 
flow produced by 

d 771 ( Po 2 



The explicit form of h follows from (5.6). 

In the above mappings through (5.8) or (5.10) the variables q 2 ,po,Pi and 
accordingly V, T are taken as constants of the motion. We can also build Hamil- 
tonian mappings which modify V or T. As a simple example, let us construct a 
transformation among the set of van der Waals equations which will make the 
law of corresponding states obvious. The critical point is given by 

V c = 36, P c = -iL T c = J^-, (5.11) 
27b 2 27b ' 

and the P, V, T equation of state depends only on the reduced variables V/V c , Pj 
P c , T/T c . We look for a mapping which leaves these reduced variables invariant. 
If moreover we wish T = —po/pi and hence T c to remain constant in time, we 
should find that the variables V = q 2 and P -1 = p\/p2 as well as the parameters 
a and b in the equation of state vary proportionally to one another. On the other 
hand, the variations of V and b entail that the contribution J? In (V — b) to the 
entropy S = q a depends on the time r. If we assume this dependence to be 
linear, in —Rt, the common dependence of V = q 2 ,b, and P -1 = p\/p2 should 
be exponential, in c~ T . The simplest Hamiltonian which achieves these goals is 

h = -Rp - P2 q 2 , (5.12) 

since it provides the equations of motion 



dp dpi dp 2 dq° dq 1 dq 2 

~T~ = ~T = °> ~T~ = P2, -r- = ~R, —r~ = 0, — — = —q , (5.13 
dr dr dr dr dr dr 



implying V cx P 1 oc e T and constant T. The left-hand side of the P, V, T 
equation of state, which is at the initial time 

{ p +yi) { y- h ^- RT= ^ ( 5 - 14 ) 

retains as expected its form in this motion of P and V, with a and b changed into 
ac~ T and bcT T . The remaining two equations of state (5.5), which involve S = q° 
and U = q 1 , also retain their form, and the choice (5.12) for the Hamiltonian 
entails that the function ^ (T) does not depend on r. Adding to the Hamiltonian 
(5.12) a homogeneous function of p and pi with degree 1 would result in a 
change of VP (T) in the free energy (5.6). 

The molar chemical potential fiN& obtained from (5.4) and (5.13) varies 
linearly as RTt along the trajectories. If we consider the two end points of the 
plateau in the isotherms which describes liquid- vapour equilibrium, the equality 
of the intensive variables T, P, /i at these points is maintained along the trajec- 
tories, so that the flow maps the saturation curve into that of the transformed 



22 



fluid. This is consistent with Maxwell's construction and with the conservation 
of areas in the P, V plane in a symplectic flow with constant T. Note finally 
that the ideal gas equation of state remains unchanged in the mock dynamics 
generated by (5.12), and that a van der Waals fluid is changed into an ideal gas 
by this dynamics in the large r limit. 



6 Metric structure 

We have not yet dealt with the concavity of the entropy S (q 1 , ■■■,q n ) as function 
of the extensive variables, which expresses the stability of equilibrium states. 
This property produces constraints on the physical manifolds M in the 2n + 2- 
dimcnsional space, that we wish to express. It entails the positivity of the 
matrix —d 2 S/dq l dq^, or equivalently the existence of a metric structure in the 
n-dimensional space q l relying on the quadratic form 

ds 2 = -d 2 S = -Y T^-dq' dq\ (6.1) 

which defines a distance between two neighbouring thermodynamic states. 

The possibility of endowing equilibrium thermodynamics with some Rieman- 
nian structure is known since a long time |l8|, but the metric (6.1) is not 
the only possible one in the framework of macroscopic thermodynamics. For 
instance, ds 2 — d 2 U also appears acceptable, although ds 2 = —d 2 S is preferable 
in the framework of thermodynamic fluctuation theory p9| . 

This arbitrariness can be lifted by resorting to statistical mechanics. The 
problem there is to find a physically meaningful distance Q between two arbi- 
trary neigbouring density operators D and D + dD. To this aim one should rely 
on the formal structure of quantum statistical mechanics. The essential fact is 
that the only quantities which behave as scalars in a change of representation 
and which have physical relevance at a givent time are (i) the expectation values 
(A) = TrDA of the Hermitean operators A in Hilbert space which represent any 
physical variable, and (ii) the von Neumann entropy 

SVn = -knTrDlnD. (6.2) 

Hence the operators D on the one hand and A or In D on the other hand can be 
regarded as belonging to two dual vector spaces. From the above considerations 
one can show [|| that the only quadratic form in dD which has physical meaning 
is, within a multiplicative constant, 



ds = —d S v n = /crTt 



dD) IdlnD 



(6.3) 



It is positive owing to the concavity of the von Neumann entropy and there- 
fore defines a natural Riemannian metric in the space of density operators D. 
The matrix elements of dD and k-QdhiD appear as the contravariant and the 
covariant coordinates of a shift. Returning to thermodynamics, we can restrict 



23 



this metric to the subset of the space D constituted by the Gibbsian density 
operators, which have a form of the type (2.1) and which represent at the mi- 
croscopic scale the local equilibrium states of thermodynamics. These density 
operators are in one-to-one correspondence with the points q , ...,q n and their 
von Neumann entropy (6.2) is identified with the entropy S {q 1 , q n ) of ther- 
modynamics. The reduction of the metric (6.3) to this n-dimensional subset 
of the space D is identical with the metric (6.1), which thus arises as the only 
natural metric inferred on equilibrium or non-equilibrium thermodynamics by 
the unique metric of quantum statistical physics. 

In this microscopic approach the volume of the system is regarded as fixed. 
It is not included in the variables q l which characterize the state of the sys- 
tem. If the set q l contains all the extensive variables, including the volume, an 
eigenvalue appears in the metric tensor (6.1), since the extensivity of S then 
implies 

y d 2 s j = 

dq i dqi q 

The distance (6.1) between two states also vanishes in the thermodynamic limit 
in a phase equilibrium situation, when these two states differ only through the 
proportion of the phases. Distances are always positive between two states 
having different intensive variables, and they thus characterize their qualitative 
differences. 

The condition ds 2 > expresses one of the Laws of thermodynamics, the 
concavity of the entropy function. In order to complete the formulation or 

thermodynamics in the spaces T, T or T, we therefore have to rewrite (6.1) 

in terms of points on a manifold M,M or M. We first note that, for a given 

manifold M, the definition (1.1) of the intensive variables r y i implies 

** = t ^ 

Thus, if we parametrize the n-dimensional manifold M associated with a given 

system in the 2n-dimcnsional configuration space T by means of the coordinates 
q 1 , q n , the metric (6.1) takes a simple form: An infinitesimal shift of the 
state can be represented either by the set of variations dq 1 , dq n which can be 
regarded as its contravariant components, or by the set —dj 1 , — dj n which 
according to Eq. (6.5), appear as its covariant ones. The distance between two 
states of any given physical system takes, for infinitesimal variations on M, the 
form 

n 

ds 2 = -J2dlidq\ (6.6) 

i=l 
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(6.4) 



which does not depend explicitly on the entropy function of the system. The 

specificity of this system is reflected by the curvature of the manifold M, itself 
resulting from the metric tensor given by (6.5). 

Considered in the whole space T, the expression (6.6), which equivalently 
reads 



2 



i (dq* + d lt ) 



2 



i=l 

defines a pseudo-Euclidean metric with signature (n, n) . It may have either sign 
between two neighbouring points of T which correspond to different physical sys- 
tems. However, the restriction of (6.6) to a thermodynamic manifold M defines 

the positive Riemannian metric (6.1). Conversely, the submanifolds M of T 
which can be thermodynamically admissible are the Lagrangean manifolds over 
which (6.6) moreover induces a positive Riemannian metric. The first condition 

ensures the existence of a function S which generates the equations of M in the 
form (1.1), and the second one ensures the concavity of this function. 

The pseudo-Euclidean metric (6.6) also holds for the 2n + 1-dimensional 

thermodynamic space T. The positivity of the metric that it induces on a n- 

dimcnsional Legendrean submanifold M of T now ensures that this manifold 
may represent a thermodynamically stable physical system. 

Similar features occur in the 2n+2-dimensional space T. The pseudo-Euclidean 
metric (6.6) can be rewritten as 

1 " 

ds 2 = — V d Pi dq\ (6.7) 

The factor 1/po ensures gauge invariance, and apart from it ds 2 is symmetric 
with respect to i = 0, 1, ...,n. The restriction of (6.7) to a n + 1-dimensional 

n 

Lagrangean manifold M where u> = ^2 Pi d<f = (which ensures the existence 

i=0 

of S and the gauge invariance) yields again the metric (6.1) on M. However, 
there are here n+ 1 contravariant coordinates dq° , dq 1 , ...,dq n for a shift in 
M, which are not independent since u = 0. The independent variables on M 
being taken, for instance, as q 1 , q n ,po, the covariant coordinates of a shift 
are dpi/po (i = 0, 1, ...,n). The distance between two points describing the same 
state, which are deduced from each other through a gauge transformation (2.3), 
vanishes. More generally two different states with the same intensive variables 

7i,—,7n ue at a vanishing distance, as obvious from (6.6) for M or M and 
from (6.7) in a gauge with constant p for M; this expresses that the metric 
characterizes the nature of the independent phases of the system, not their 
sizes. 
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Altogether the physical manifolds M in T should not only satisfy the condi- 
tion lo — 0, which ensures that the thermodynamic identities are satisfied, but 
should also be such that the bilinear differential form (6.7) induces on them a 
non-negative Riemannian metric ds 2 , which ensures thermodynamic stability. 

Accordingly canonical transformations 7i which map M onto M' in T pre- 
serve thermodynamic stability only if they satisfy some inequalities. Using Eqs. 
(3.8), we can compare the metrics at corresponding points through 

) E [ejw*>* d « J + WW dq dq ) ' 
^ E ( irnn d ti dq3 + irnrrM dp A ■ ^ 

The mapping should be such that, if the points q l ,Pi and </\p- (i = 0, 1, ...,n) 
and their variations lie on the manifolds M and M', respectively, both ds 2 and 
ds' 2 are positive. In a continuous transformation generated by h, the metric 
evolves according to 

d (J 2\ 1 dh j2 1 / d 2 h d 2 h , tJ A 

* {ds) = P~oW ds + 7, £ n (^-^ ^ - w d9 d(Z J ■ (6 ' 9) 

A control of the sign of the metric through Eqs. (6.8) or (6.9) sets up conditions 
on TC or h. The occurrence of a zero eigenvalue in the metric tensor during the 
evolution in terms of the mock time r indicates the appearance of a critical 
point in the equations of state thus generated. 

7 Conclusion 

The extension of the n-dimensional configuration space q 1 , ...,q n associated with 

a thermodynamic system, successively to the 2n-dimensional space T including 

the intensive variables 7i , ■ ■ ■ , 7 n i tnen to the 2n+ 1-dimensional space T including 
in addition the entropy q° , finally to the 2n + 2-dimensional space T where the 
7j's are replaced by the variables po,Pi, ■■■,Pn including a gauge arbitrariness, 
presents practical advantages. It is adapted to changes of variables among the 
various physical quantities. It sets into the same framework all the systems 
which have the same degrees of freedom q 1 1 ...,q n but different equations of 
state. Moreover it discloses a rich mathematical structure that characterizes 

the geometry of the spaces T, T or T, and of the manifolds M, M or M which 
represent a thermodynamic system. 

The existence of thermodynamic identities is reflected in the Legendre struc- 

ture of any M, characterized by the identity u) = dq° — li^q 1 = 0, or equiv- 

i=l 

alently in the Lagrange structure of any M or M characterized by the more 



ds 2 = 
ds' 2 = 
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symmetric identity u = V ; A<£ = 0. The conservation of ui in transforma- 

i=0 

tions which map the manifolds M onto one another induces for T a symplectic 

n 

structure based on the 2-form du; = dpi A dq l , where the variables q l and pi 

i=0 

appear as conjugate in the sense of analytical mechanics. Such mappings are 
then canonical, and can be generated by Hamiltonian flows. The quantities q 1 
play the role of position variables, the quantities pi or -f i that of momentum 
variables. The Hamiltonians involved should be homogeneous with degree 1 in 
the variables pi so as to preserve the gauge invariance of the manifolds M. The 
mappings such that the trajectory issued from any point of some given manifold 
Mo lies in Mo are those for which the Hamiltonian vanishes over Mo. Finally 
the stability of thermodynamic states is reflected in the fact that the bilinear 

n n — 

forms — ( h i dq l or p^ 1 dpi dq l should induce on M, M or M a Riemannian 

i=l i=0 

structure with a non-negative metric. 

We have shown how Hamiltonians can be constructed, either for real motions 
in thermodynamic processes, even when they are dissipative, or in fictitious 
dynamics which consistently generate equations of state for different systems. 

The above considerations should have more than a formal value. We have 
worked out a few elementary examples for the sole purpose of illustration. More 
realistic problems should benefit from the existence of an ordinary Hamiltonian 
structure, which for instance readily generates variational techniques, and which 
anyhow is formally simpler than the currently used contact structure. 

We wish to thank the Total Fina Elf group for support and permission to publish. 
We are also grateful to the referee for his thorough screening of the manuscript. r 
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